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ABSTRACT 


The vibration analysis of various types of 
plates namely rectangular (with and without cut outs), 
annular, annular sector, circular and rotating discs has 
been done. For the case of rectangular plate, a four 
noded rectangular element with foxir degrees of freedom 
at each node, namely lateral displacement, two bending 
rotations about two co-ordinate axes respectively and 
twist has been used, 

A circular sector element with same number of 
degrees of freedom and an annular element ( semi analytical 
method) have been used in case of circular plates. The 
annular element takes inpD consideration the thickness 
variation in tl^e radial direction. The vibration analysis 
of rotating discs has been done using the semi analytical 
method. 

Nowadays the higher stiffness to weight ratio 
of orthotropic plates is drawing much attention. Hence 
throughout the formulation and analysis of results , 
orthotropic plates are given enphasis. 



CHAPTER - I 


INTRODUCTION 

1..1 GENERAL- 

One of the cornmonl-y used structural conponent 
in the industrialised vjorld is a plate. In many of the 
design problem specifications merely ensurirg that plates 
will v;ithstand the applied static loads will prove to be 
inadequate since a variety of structures used in land/ 
sea, air and space are subjected to dynamic stresses 
and displacements induced by periodic forces acting on 
the lateral surface of the plate. Random forces are 
expected, for example, on the surfaces of the plates 
exposed to tangential gas flow as in aircraft cjomponents 
or stationary structur-es exposed to high wind velocity. 
Such forces can also be expected when a fluid flov7 takes 
place along the plates as in the case of rectangular 
plates used as the hull of a ship. 

Periodic forces are likely to be experienced 
when plates form a part of the rotating machinery operat~ 
ing at a specified speed. To ensure mechanical integrity 
it is also necessary to predict the influence of stresses 
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in the discs which are due to the centrifugal force 
arising in the rotating disc due to its own mass. 
Turbomachines and superchargers are typical examples 
where such rotating discs are used. 

Usually little can be done to change the nature 
of the driving forces. Kence the most commonly used 
design technique for solving such problems where matching 
of the driving forces and. natural frequencies is to be 
done is to alter the plate geometrq?' or boundary condi- 
tions/ so that the natural frequency zone is away from 
the frequency zone of the driving energy and hence the 
resonance is avoided. The importance» of the knowledge 
of the natural frequencies of plates in such applications 
is obvious. 

1.2 LITERATURE SURVEY 

In literature one finds an extensive treatment 
of the problems on the flexural vibration of isotropic 
plates - both rectangular and circular. The piublem of 
vibration of isotropic plates having free edges was 
investigated long back by Kirchhoff, Lamb and Rayleigh 
using Poisson-Kirchhoff theory. Timshenko [l] used 
energy method for solving the case of isotropic plate 
with clamped edges, Thein Wah [ 2 ] investigated the 
case of simply supported edge using Poisson-Kirchhoff 
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theory for the basic foirmulation. Numerical results are 
also given for the frequency parameter, for modes of 
vibration consisting of nodal circles and diameters. 

An extensive study of uniform annular isotropic plates 
has been done by Vogel et al. [3]. Here the classical 
theory of flexural motions of elastic plates was used to 
determine the natural frequencies for various boundary 
conditions of annular plates. An investigation was also 
m.ade into tlie effects of Poisson's ratio on frequency 
parameters by talcing different values of Poisson's ratio. 
Its effect has been found to be insignificant on natural 
frequencies. Hence in reference [3] the analysis has 
been done taking the value of Poisson's ratio as 0 . 3 . 

Blevins [4] has given numerical results for 
rectangular , annular and circular plates for various 
boundary conditions. Zinckiewic? [s] has given applica- 
tions of finite element method to the vibration analysis 
of plates. The common elements used are triangular and 
rectangular [d] ones. Triangular elements are more use- 
ful since they can be used to represent curved boundaries. 

Ergutoudis [?] has developed quadrilateral 
element v/ith curved sides. Olson et al, [s] have deve- 
loped two plate bending finite elements in polar co- 
ordinates to study the vibration of circular plates. 
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The first element has three nodal corners and is in the 
shape of the sector of a circle. The second element has 
four nodal corners and foirms an annular sector. The 
transverse displacement and two rotations about the two 
co-ordinate axes respectively are the three d€ 3 grees of 
freedom used at each node in both the cases. Non-dimen- 
sional frequency parameter for a circular plate with 
clamped/ free and simply supported boundary conditions 
are tabulated and compared with the exact values. Kirhhope 
et al . [9] have used semi analytical method to study the 
dynamic behaviour of isotropic circular plates and rotat- 
ing discs for various boundary conditions, Kultar Singh 
[10] and Benjamin [ll] have also used the annular element 
to study the dynamic behaviour of plates . 

Even though the literature contains various 
types of solution techniq^ios of vibration analysis of 
isotropic plates /that of orthotropic plates needs some 
more attention, Kirmser et, al. [12]/ Pandalai et al. 
[13]/ Minkarah et al, [l 4 ]/ Joung [isj, Bellini [is] 

Huang [l?] and Nowinski et al, [is] have attempted some 
problems on orthotropic bodies subjected to dynamic 
loading. The existance of material singularity has been 
reported by many of thai;!. This is duo to fact that 
circumferential and radial rrwdulai can not be different 
at the centre. In order to overcome the above problem 
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some have assumed an isotropic core of small radius at 
the centre. 

Ramaiah and Vijayakumar [l9] have found out the 
lox-irest eigen values for various houndaary conditions using 
Raylsigh-Ritz method x^ith simple polynomials as admissibl 
functions. In an atterapt to find out the natural fre- 
quencies of higher modes of vibration / Ramaiah et al . 
[20] have used the classical Rayleigh-Ritz method by 
introducing the co-ordinate transformation. They have 
obtained the natural frequencies corresponding to the 
asymmetric modes (modes X'^ith one nodal circle and tx-jo 
nodal diameters). Hox^/ever application of tliis method 
is generally cumbersome since the method involves solu- 
tion of a large set of simultaneous equations for higher 
accuracy’' of results. They have given the results for 
free-free/ clamped-free and simply supported - free 
annular plates of radii ratio 0,5 and 0,1, Regarding the 
vibration of rotating discs Evessham's [21I work can be 
mentioned here where a series solution has been obtained. 

1.3 PRESENT vroRK 

The study of nature of the natural frequencies 
of plates (mainly orthotropic) for various boundary con- 
ditions with vaarying thickness and rotating discs is 
the aim of the present work. Using rectangular and 



6 


annular sector elements (described in Chapter II) vibra- 
tion analysis of plates ’ ith geometrical nonlinearity 
like plates with cut-outs can be m^de. It can also take 
into account different t}rpes of loading. 

For plates where there is linearity in material 
property as well as symmetry in geometry semi analytical 
method is more efficient reducing the problem to a one 
dimensional one. In the present vrork vibration analysis 
of orthotropic plates with varying thickness and rotat- 
ing orthotropic discs has been done using the above 
method (Chapter III and IV). For uniform thickness 
plates good results are obtained v/ith few elements ( 2 
to 4). But slightly more number of elements (8 to 10) 
are required for plates v/ith varying thickness. 



CHAPTER - II 


VIBRATION /^JALYSIS OP RECTANGULAR AND ANITULAR PLATES 


2 . 1 INTRODUCTION 

For vibration analysis of rectangular plates/ 
recta33gular elements have been used. A typical rectan- 
gular element is shown in Fig. 1. The element chosen is 
a conforming element with four degrees of freedom per 

node, they being deflection w, two rotations W , W and 

X y 

twist This rectangular element with these sixteen 

degrees of freedom is compatible in displacement and its 
first derivatives i.e. it is a Cl element. For circular 
plates an annular sector elsnent, shown in Fig. 2 has 
been chosen with same sixteen nodal degrees of freedom 
in polar co-ordinates. 

2.2 ASSUMPTIONS 

The present formulation is based on the thin 
plate theory (for small deflections) which mabes the 
following assumptions [ 22 ]. 

1, The plates deform through flexural deformation only. 
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2. Deformations are small in compariOTr. with the plate 
thickness. 

3. Normals to the mid-surface of the undeformed plate 
ronain straight and normal to the middle plane during 
bending . 

4. Rotary inertia and shear deformation are negligible. 

5. 'Hie inplane stresses (mombrane stresses) are negligible. 

2.3 formulation OF THE PROBLEM 

Following formulation is easily available in 
literature and is given here for easy reference. For an 
elastic linear material strain ene3?gy is 


rHfcM 

fl 

1 [ ej 

V 

[d' ] i e [ 

L J 

dv 

(2.1) 

and the kinetic 

energy is 




T = 1 

/ p h ^ 

A 

dA 


(2.2) 

(e) 

Assuming W = 

LnJ ( w 
L 

} (ne) 


(2.3a) 

Using eqn, (2.3a) J = 

= [b*] 5wt 

(ne) 

(2.3b) 


= 1 [wjCi'e) j [B'f [d'] [b'] dv Jw 

_ J 


(2.4) 

= -J / P h J n} [n] dA ) w 

.... (2.5) 
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Now applying the Hamilton principle 

(T ~ U) dt = 0 
t^ 


( 2 . 6 ) 


Substituting eqns. (2.4) and (2.5) in eqn. (2.6) one gets 


^ W* + [k]^®^ I W [ 


t ^ 


n°i - 

(2.7a) 


where [m]^®^ = / p h ^ N c [Nj dA and 

A 


1 J 


= / [B'f [d'] [b'] dv 


(2.7b) 


V 


Asstuning harronic motion, eqn. (2.7) becomes 

( [k]^®^ - [m]^®^) J W V = -<0 \ (2.8) 

L 

which are the equations of a motion for any elsnent. 


2.3.1 Rectangular Plate 

For a rectangular plate s-^rain matrix [s] is 

[22] 


le J = L® ® 7 J 

L J x Y xy-* 

= I - z,^ -z.-i 

♦vx OY 


. 2Z 

d^dY 


(2.9) 

Using eqn. (2.3) and eqn. (2.9), [s'] matrix for rectan- 
gular elements becomes 
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[b'] = Z 




^XX' 


- LH,yyJ 
-2LK,^J 


= z [b] (2.10) 


i(e) 


Using eqn, (2,10), stiffness matrix ^2. 7b) 

becomes 


[k]<®>=// [bF [d] [b] dx ay 
For orthotropic plates [d] matrix is 


[d] 


X 

D. 


D. 


D 


0 


O 

0 


D 


xy 




""x = 12(1-V^ Vy ) 
E 

y 12 (i-v^ V ) 


=: ^X ^y ' 
Eqn. (2.12) reduces to 


D 


_ G h /12 
xy ” 


[d] = 


E h' 


1 

V 

o 


V 

1 

0 


0 

0 


1-v . 


(2.11) 


( 2 . 12 ) 


(2.13) 


12(1-V^) 

r i(ef ^ ^-n (2.7b) becomes 

Consistent mass matrix of eqn. v 

l(e) _ rt 0 T, \ KT : (Nj dx dy (2.14) 


= // P h < N I 

K jf 



2.3.2 Shape Functions 
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As said plate bending element needs the compa- 
tibility Of displacements and slopes. The shape func- 
tions satisfying these conditions can be obtained from 
the follov;ing matrices product [23], 


( e) 

^ = IN, N 

lx 


N 


3x 


N. J 

4x 


W 3 

W 

^3 

^2 

w 

^2 

^x 

"^3 

W 

xy3 

^2 

W 

xy^ 

^4 

Wy 

^4 


W 

^1 

^x 

^4 

w 

^^4 

W 

w 

xy^ 



where and are the standard beam shape functions 
and are given in Appendix-A. Rearranging eqn. (2.15) 


= LNj N2....NJ5J ; 


= LNJ 

(2.16) 


, , 

where |^W| is nodal displacement matrix (w, 

' ^xy. ' i = 1#4) 


'1 ' \ ^ 


Yi 

It may be noted that shape functions [Nj are represented 
interms of local co-ordinates ^ and usirg- the 
transformation 


^ = 


3C - X 


X 


T1 = 


y - y 


( 2 . 17 ) 
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2.3.3 Annular Plates 

For circular plates the strain matrix lS-J 
from ref. [ 4 ] is 


L^J = I -2 


■dfw 

dr2 


-ZCl'iW + _2z(i 


r ;)r ^ ^ , 2 


1 &?w 


- 1 '' 1 - 2)1 

^orO© _2.)© 


r 

... ( 2 . 18 ) 

said earlier sector element shown in Pig. 2 has been 
chosen for circular plates. Requirements of shape funo- 
tions being same a= oarller case,the shape functions used 
for rectangular plates can be used. Usi.^ eqns. (2.16) 
and (2. 18), one gets from eqn, (2.3b) 


[bO = 2 


TT 


- + ^2 K,ggJ 


-2f~ N, 

"r 'r© 


1^ 


= 2 [b] 


2 


( 2 . 1 ?) 

Using eqn. (2.19), stiffness matrix of eqn.(2.71o) 

becomes 


= s: [sf [d] [b] r dr ae 
where [d] matrix for orthotropic plates is 


( 2 . 20 ) 


~ 12(l-v^ Vq) 


E 


E V- 
r © 


\ ^© ^ 


0 


0 


0 

0 


G^Cl-V^V^) 


• •. m 


( 2 . 21 ) 
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and for isotrop-ic plates it becomes 


[d] = 


E h‘ 


1 

V 

•i 

i 

I 0 


V 0 




1 0 
0 


1-v 


( 2 . 22 ) 


12 ( 1-v ) 

Consistent mass matrix of eqn. {2.7b) becomes 


= P h // [ r dr d© (2.23) 

In order to do the numerical integration, the sector 
element has to be mapped into a square of side 2 units. 
This is done using the following relations [25] 


5 


(r - b) 


1 and n = 


2 (© - 0 ^) 

-T57-^) 


(2.24) 


From eqn. (2.24)/ r can be written as 

r = I (1 (a - b) + b (2.25) 

2.3.4 Numerical Integration 


For this Gauss Legendre quadrature has been 

used. 

Rectangular Eleme nts; For evaluating stiffness 

and consistent mass matrix one needs 

the second derivatives of the displacement with 
reference to global co-ordinates and global 
differential area in tearms of local co-ordinates 
i: , n / This is done as follows’ 
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First, derivates of the displacement with 
■eference to global co-ordinates are given by [s] 


(2,26) 


where [j] is the Jacobi matrix and is given by 


f a-i 


afi 1 

-:■ ^ 
j ANi i 

(f . 

1 — 1 

1 j 

/■ dx ^ 
'■ ^ f 

: ^n 


j dy 

L vJ 

1 






[j]= 


(2.27) 


'■^ 2 ^. w£v: 
jrn an 

Differentiating eqn. (2.26), (first eqn. w.r.t. t 


and second eqn. w.r.t. 7^ ), 

■ j ^ 

X ^ ■ X (.jx 


*% U ; 

lift 

, on^ ! 


I 3 y] 

[j] J dx^ o e ^tJ 

Vo " o ^ 


dx^ ay^ 

1 

I I 

_ _ I >. _ 2 . 


= [a] [.] ^ 

O Nf 


(2.23) 


Differentiating first eqn. of eqns. (2.26) w.r.t. 
Tf' / one gets 

A^N. „ ^ '^2 jN. 

tf i. _ (■ u_x itjr 4X> 11 4. Js'LJL. ^L-i. 

A? 5n q>i' An dn dx^y ao^''-n'^ 


o^i 


cj ^ AT) 
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For rectangular elenents this equation siinplifies 
to 








(2.29) 


al bi 

For this rectangular element, using eqn. (2.17), 
[j] matrix given by eqn. (2.27) becomes 


[J] 


0 


0 


L : 

Vi 


(2.30) 


It is well known that 


/ / dx dy = 


J 


IS 




an 


(2.31) 


Sect o r El emen t s; For these elements the above deri- 
vation is valid except that x and y should be 
replaced by r and 0- ; and Jacobian matrix 


becomes 


V] = 


and 


1 


a— b 
2 ' 


0 


® 2 ~® 1 ! 




// dx dy 


;/ r dt dn 

7 


(2.32) 


(2.33) 


2,4 SOLUTION PROCEDURE 

Stiffness and consistent mass [m]^® 

matrices have been calculated usir^ the numerical 
integration explained above. It may be ixjted that 



the elements used are subparametric elements. These 
elemental matrices are assembled as banded matrices 
after applying the relevant boundary conditions. 
Eigenvalues are calculated from these using the NAG 
subroutine available on DEC 1090. 



CHAPTE R I II 

analyst S_ OF CIRCULAR AND ANNULAR 

fpKpMC v^ying~''to 

3.1 INTRODUCTION 

The standard finite element method described 
earlier is no doubt capable of dealing with any plate 
problem. Nevertheless the cost of solution increases 
rapidly with higher number of elements which are required 
for higher accuracy of the results. So in the case of 
circular or annular plates where alonq ©-direction there 
is no change either in geometry or material property 
semi analytical method can be used [5,9 ]. It reduces the 
two dimensional problem to a one dimensional one and hence 
very economical. Here we express the field variable in 
G“d.irection in Fourier series, 

3.2 FINITE ELEMENT FORMULATION 

Strain energy of a finite element is given by 

i / [e] [d'’] |e I dv (3.1) 

V • V J 

(e) a-T- th 

Taking the displacement W" corresponding to the m 

harmonic over the element (Pig. 3) as 
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or 


where 


and 


rCe) 


^^l+9-2^+^3^^+a-4r^) oos m © 

= LSj [r] cos m © ) W f 
Lsj r: [i r r^ r^J 


^ W , 
i. J 


Lwj 


= LW^ ^2 J 


1 “ "2 
Cr] is given in Appendix-B. 


Using eqn. (3.2), strain matrix becomes 


( 3.2) 


or 


(' ) 

(_z) 


r 


cos m © 


m 


L? ='r - — 


sj cos m © 


2L“ “ S,^ + Sj sin m © 
r 


t ) 

) et = (-Z)i 

^ J 


[r] 


2 cosrn© 6 r cosifi© 


rn ^ 1 m'^ 

~ --gcjosm© ) cosm© (2-m^)cosm0 r( 3-in^) cosrn© 


I “slnn:© 

r 


0 


-m sinm© 


-2m r sinm© 


or 


X [r] i MV 

L J 

,^ej= (-Z) [pO [r] 


(3.3) 

SUbstitutir^ egn. (3.3) in eqn. (3.1)/ and integrating 
’w.r.t. Z © one gets 

(e) 1 ,i..i(ne) r.,iT ^-^3 r.,iT r.,1 , r..-! v„! (ne) 

;{_W5 

(3.4) 


= I [r]^ / [p]"^ [d] [p] r dr [r] [w { 
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where 


TaTPV-Vg) 


E V_ 

r 0 


V 


G _ { 1-v v_ 
r0 r 0 


(3,5) 


C =r 2^ 


C = T 


when m - 0 


when m > 1 


(3.6) 


( 2-m^ ) ( 3-m^3 


(3.7) 


Eqn. (3*4) gives us the sUff ness matrix 


= [Rf j h3 [rf [d] [p] r ar [r] 

h 


= [nf [r] (3.8) 

Coefficients of are determined by closed form 

integration and are given in Appendix-B, 

Kinetic energy of a finite element is given 


= - IS 9h. r dr d0 

Substituting eqn, (3.2) in eqn. (3.9) one gets 


(3.9) 
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^ 1 C ^ I 

from Whidi one gets the cx>nsistent mass matrix [M]^®^as 
[m]^®^ = c [r]*^ /Ph LsJ r dr [r] 

= g [Rf [r] (3.10) 

Coefficients of [m]^® are given in appendix~B, 

3.3 SOLUTION TECHNIQUE 

Elemental stiffness [k]^®^ and consistent mass 
[m 3^®^ matrices are assembled into global matrices. 
After applying the boundary conditions, the eigenvalues 
are obtained using the IMSL subroutine EIGZP. 



CHAPTER - IV 


FREE VIBRATION OF R OTATING _pi SCS 
4.1 INTRODUCTION 

The vibration of rotating discs is of nrudh 
practical significance. In cases* 1‘ilce symmetric discs 
■used in axial flow turbomachines / formulation of the 
previous chapter can not be used. Now centrifxigal 
force becomes an important factor. Its contribution 
to stiffness matrix must be includejd.To do so one needs 
the inplane stresses. First these inplane stresses 
arising due to the centrifugal force are determined. 

In the present work ther 3 stresses have been calculated 
using finite element method. The formulation is appli- 
cable to both isotropic and orthotropic discs. These 
inplane stresses make the disc more stiffer to flexural 
vibrations and hence natural frequency increases. Add- 
ing the additional stiffness matrijg obtained due to 
these stresS/to the original stiffness matrix given in 
Chapter-Ill we get the stiffness matrix for rotating 


d i SCS 
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The formulat: jn below is a ^ jneralised one 
which can be ’.ised fox* variable bl '.ckness discs (isotro- 
pic and ortho tropic) . 


4 . 2 P0RI4ULAT10N 


The strain energy due to the inplane stresses 
during flexural vibration [ 9 ] is 


U 


(e) 


P 


1 2v a 

~ S J 

"oh 




or 


_9 

.2 


/ U \ h 2 -} 
j 


h r dr d0 

(4,1^ 


(e) 

Substitutlrig W frxjm e'ln, (3.2)/ we c et stiffness 
matrix as 


= c [af I h fa-r is,p+ffg ^ (s|lsj)r teW 

= c [af LS] f.i.2) 

where C =: 2v for m ~ O ( ■■ 3) 

C = V for m > 1 

The coefficients of the matrix [h are calculated 

by closed form integration and are given in Appendix-C. 


4.3 DETER14INATI0N OF liiPLANE STRESSES 

The inplane stresses and needed for 

eqn, (4.2) can be obtained using the same annular 

(e) 

element. The strain eneirgy U) of the annular ele- 


ment is 
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- -2 2"^ / h L^J [d] r iir 


(4.4) 


where 


;e 


[d] = 


1-v 

r O 


r €H 


E 

r G © 


I and L^' J = 


, du u j 
‘■dr r-^ 


( e) 

Taking u" over the element as 


(4.5) 


u = Ca^ + a2r + + a^r“) 

= LHJ |u? 
where [nJ = l_Sj [r] 

Substituting eqn. (4.6) in ean. (4.4) one gets 

u{®> = 2^ Luj'“> / h [sf [d] [b] r dr 

r 1 ... (4.7) 

i Ln, j ! 

where [b] = '> (4.8) 

I- J 

Eqn. (4.7) gives the stiffness matrix 




u 


= LN, N 3 NJJ. - 1 , 

' i 

U ! 


(4.6) 


= 2m / h [b]"^ [d] [b] r dr 

= [Rf (4.9) 

Coefficients of ma-fcrix are obtained using closed 

form integration and are given in Appendix-C. 
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The load matrix due to body force (centrifugal 

force) is 

4 = 2-^ph I tNj dr (4.10) 

where P is the mass density, U is the angular velocity''. 
These elemental matrices are assembled and displacements 
are obtained. Using these displacements, stresses 
and O ' 0 are .calculated using eqn. ( 4 . 5 ). 

Usir^ these values of ^ j- and , stiffness 
matrix for the disc is calculated. VJhile cal- 
culating the [k stress distribution inside the 

ir 

element is taken linear. 

i .e. 


where 


Or = El + r 



By adding [k to of plate ( Chapter-Ill ) one 

Ir 

gets the stiffness matrix for the rotating disc. This 
is solved for finding out the natural frequencies of 


the rotating discs. 


Results for isotropic plates are also calcu- 


lated usir^ the closed form solution of stresses of 
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CHAPTER - V 


fiESULTS_ AND JDI _SaJSSiqN 

5.1 RECTANGULAR ELEI'-IENT 

The finite element program developed has been 
verified for convergence and accuracy. Natural freguend 
of square plates with all edges free and clamped are 
determined and compared with results from ref. [ 4 ]. The 
convercence is checked by varying the n-umber of elements 
taben (Table la). The natural frequency of a plate with 
a centrally located squre hole is also determined. With 
12 elements the answer is found very close to the results 
from ref. [ 4 ]. (Table lb). In all the cases, the natural 
frequency is expressed in non-dimenoional form 

a = X 1^ ( P h/D) 

o 2 

where D = E hVl2(l-v ) and 1 is the width of the plate. 

5.2 ANNULAR SECTOR ELEMENT 

Here the elejment is used for the determina- 
tion of natural frequencies of ahannular sector plate 
with all edges clamped having a radii of 0.4 an<J s-ecto- 
fial angle of 90°., The corwergeiMie can be verified by 



29 


Comparison of the frequency parameter Tv (Ph/D) for 
a rectangular plate with answers from ref. 





grid size 

B-C 

Mode 

2x2 

3x3 4x4 ^ 


1 

7.54 

13.67 

13.49 

P"F 

2 

20.23 

19.32 

19.79 


3 

43.33 

24.76 

24.43 


4 

56-09 

35.43 

35.02 


1 

36 .75 

36.05 

35.99 

2 

C-C 

76.01 

'74.03 

73.41 

3 

76.33 

74.03 

73.41 

4 

113.61 

109.43 

103-30 


Table la. 


Comparison of frequency parameter for a clamped square 
plate with a square hole of -side lergth 1/2 located at 
the centre with answers from ref. [47*" 


B-rC grid size F.E.M. jf 

C"C 4x4 66.09 62.40 


Table lb 


Comparison of the frequency parameter K R^(ph/D) of a 
clamped annular sector plate with answers from ref. [26?^ 


Mode 

- I-llSifiA 

4x4 

size 

5x5 




. .. 



1 

5 3.63 

5 3.41 

52.70 

2 

91.90 

91.60 

33.30 

■ 3 

123.10 

122.07 

125.10 

4 

153.70 

152.76 

140.00 

5 

162.30 

157.12 

170.40 


Table Ic. 


2 1 '^2 

Comparison of the frequency parameter X R 2 (ph/D) ' for 
F~F annular plate v/ith answers from ref. [ 4 ]^ 

Mode T~x'2 X- 

1 _ 4.51 4.23 

2 12.10 9,32 


Table Id 
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comparing the results tatalated with 4x4 grid and 5x5 
grio (Table Ic) The answers from are expected 

to be higher than the exact values. Here it is violated 
tor some of the modes. This may be due to the unegual 
size of the elements used. By increasing the number of 
elements /required accuracy in answers can be obtained. 
Natural frequencies of a free— free annular plate is 
tabulated in Table Id, 

5 . 3 ANNULAR ELEMENT 
5.3a Isotropic Plates 

Because of the non-availability of the results 
for orthotropic plates for all the boundary conditions, 
first the results obtained by the semianalytical method 
for isotropic plates ara compared with exact values 
wherever available as a verification of uie program. In 
other cases only the convergence of the results is verified. 

The vibration analysis of a solid circular plate 
can also be carried out using this method mabing an assump- 
tion that a small hole (Rj^/^^2 ~ 0»001) exists at the 
centre. It is found that the natural frequency of an 
annular plate with radii ratio 0.01 itself is very close 
to that of a solid plate. As the ratio decreases further 
the frequency parameter remains almost stationary . These 
results are tabulated in Table 2. 
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Frequency parameter X ( f^h/D) for a solid circular 


isotropic plate 

for 

various 

boundary 

conditions with 

answers from ref 

p 

1 

• L_ 






Ne = 

4 





./_^2 


B-C n 

m 

0.01 

0.001 


1 

0 

9 .00 

9.00 

9.03 

1 

Free 

1 

20 . 43 

20 . 43 

20.52 

0 

2 

5-35 

5 . 35 

5.25 

0 

3 

12.43 

12.44 

12.23 

0 

0 

4.93 

4.93 

4.97 

• 0 

1 

13.90 

13.90 

13.94 

Simple Q 

2 

25,62 

25.62 

25 . 65 

1 

1 

29.74 

29.74 

29.76 


0 

10.21 

10.21 

10.22 



21.26 

21.26 

21.26 

Clamped q 

2 

34.91 

34.91 

34.33 

1 

1 

40.24 

40.24 

39.77 


Table 2 
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In case of anni'-lar plates the variation of 
-crecruency paratmter with radii ratio is shown in Figs. 

4 and 5. In case of a free— free plate the freguencv 
decreases with increase in radii ratio. In all otlier 
cases considered it goes on increasing with the radii 
ratio . 

5.3b Orthotropic Plates 

The orthotrxjpic plates are studied with three 
cases of boundary conditions namely simple-clamped, 
clamped-clamped and free-clamped. The above combination 
is chosen because of the availability of i±ie results for 
the above mentioned boundary conditions [ 20 ]. The results 
are tabulated in Table 3. 

OrthotTOpic plates having different ratios 

are studied keeping constant. Prom Pigs. 6, 7 and 8 

9 

it is clear that the fresq'^^e^cy parameter a = X (R^^) 

1 ^9 

( t'h/D^) increases with D„/D„. But variation is very 
small. This shows tliat the material property has the 
maximum influence on a compared to and It can 

also be observed from Figs. 6, 7 and 3 that a increases 
with ratio linearly. 

The natural frequency of annular plates (both 
isotropic and ortho tropic) with varying thickness is 
also determined using the same element. Here a linear 








Frequency parameter 
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Comparison of the frecruenc^r parameter X R^Cph/D ) with 
the answers from ref. [ 20 ]^ 




-Vhr. 

2.0 

D /'D 
re r 

= 0 . 35 V g 

= 0.3 








B-C 

n 

m 

1 ^ 
CD 

0.1 

0. 3 

0.5 




1 

6,516 

6.220 

5.570 


0 

2 

2 

6.502 

6.210 

5.567 




4 

6.497 

6.205 

5 .666 




1 

11.599 

11.383 

13.453 

F-F 

2 

0 

2 

11.533 

11.337 

13.444 




4 

11,513 

11.303 

13.436 





11.470 


13.503 




1 

15 .365 

15.667 

15.0 36 


1 

3 

2 

15.675 

15 . 620 

15.025 




4 

15.672 

16.610 

15.021 





15.760 


15.600 


1 6.261 6.516 7.294 


0 

0 

2 

6.213 

6.459 

7.269 



4 

6.20 3 

6.449 

7 .226 




6.250 


7.300 



1 

15 . 43 3 

14.774 

13.790 

S-F 0 

1 

2 

14. '^23 

14.423 

13.663 



4 

14.900 

14.330 

13.654 



1 

31.694 

23.927 

26.602 

0 

2 

2 

23.803 

27.740 

26.123 



4 

28.650 

27.699 

26,033 




28 ,500 


26.210 



1 

11.826 

13.252 

19.156 

0 

0 

2 

11.700 

13.136 

19 .039 



4 

11.663 

13.107 

19.075 




11.690 


19.075 


1 

24.142 

22 . 3*0 

24.205 

0 

1 

0 

1 

2 

22.436 

21.656 

23.954 



4 

22.353 

21.537 

2 3.920 

5f 


1 

5 3.137 

40.619 

37.205 

0 

2 

2 

39.613 

37.103 

36.159 


4 

33.394 

36.954 

36.045 




33.300 


36.090 


Table 3 



Frequency parameter 





Radii ratio (R^/R2) 

7 Variation of the frequency parameter wit 
fodis ratio for differont Dq/Dr ratio for a 

C"C plate : : 









A gA 
^ ■ 


W 


^'•4 


y 


Rcdii ratio 


Q -5 


!:ion of tho frequency pcromeier wi 
ratio for different (De /Or) raises (or 


4C 


variation in thickness ong the radial direction is 
considered, 'fhe thickness at outer radius R., is kept 
constant (n^sO.l cm). The thickness at the inner radius 

varied from 0,03 to 0,12 cm in steps of 0.01 cm. 
In this case fraouency parameter cc is '.'xpressed as 
ot. = R^ (l2t ( l~v^ VQ)/(h 2 Pigs. 9 and 10 show 

the . variation of a v;ith radii ratio for different thick- 
ness ratios. Prom Pig, 11 it can be noticed that varia- 
tion with thickness ratio is almost linear when thickness 
variation is small. 

I'Jhen and ~ 1 it is the case 

of isotropic uniform thickness plates. Hence in this, case 
values are compared with exact values (Table 5). It is 
a verification of the program for its dependibility. In 
Table 4 and 6 convergence of the results for vari^ing 
thickness plates is verified. 

Table 1 , 3 and 9 give the values of fjrequency 
parameter for different and ratios for simple- 

clamped, clamped-clamped and free-clamped plates respec- 
tively. 

5.3c Rotating Isotropic Discs 

The dynamic behaviour of rotating discs are 
fhe two most commonly used boundajry condi- 
tions. They are inner edge clamped or free and outer 



pc?r Qi 





f- 


A 

<• «. 


I. 


h2 


frequency 
p-Gte with 


pcrcrn’t'i;-' 
th:ckri^..^s rcliD C-B 








44 


Verification of the convergence of frequency parameter 

\ for orthotropic plates 

with linearly varying thickness 

~ ^rs'^^r “ =0-3 

^ 1^2 = 







■ “rTTr^ 




B-C 

n 

in 

Ne 


^ 

X z 





0.1 

6.2 

6.3 

6.4 

0,5 


0 

0 

3 

10 

3.74 

6.22 

3.39 

11.77 

17.43 



3,76 

6.22 

3.39 

11.77 

17.43 


0 

1 

a 

4.45 

6.40 

3.79 

12.34 

18.14 




10 

4.47 

6.40 

3.79 

22.34 

13 . 14 

F~C 

0 

2 

8 

6.66 

3.34 

10.77 

14.46 

20.46 




10 

6.66 

3.34 

10.77 

14.46 

20.46 


0 

3 

8 

12.83 

13.37 

15.56 

18 .89 

24.77 




10 

12.87 

13.72 

15,56 

18.39 

24.77 


0 

0 

8 

32.95 

47.49 

62.71 

36.05 

124.25 




10 

33.01 

47.51 

62.71 

86.21 

124.18 


o 

i 

8 

39.13 

49.80 

64.61 

37.37 

125.64 

r-r 


j. 

10 

39.21 

49.78 

64.61 

87.3 3 

125.57 


0 

9 

8 

43.49 

56.34 

70.42 

92.57 

129.33 



10 

43.48 

56.82 

70.42 

92.74 

129.75 


0 


8 

63.41 

68.39 

80.29 

101.05 

136.37 


3 

10 

6.3.43 

68.86 

80.29 

100.86 

136.79 






^l'^2 

=-0.9 




0 

0 

, g.,. 

10 

3.56 

3.56 

~T,2T 

5.22 

~7.03~ 

7.03 

9 .3“6 
9.86 

14.60 ' 

14.60 


0 


8 

3.78 

5.36 

7.36 

10.34 

15 .20 

1? r' 

1 

10 

3.79 

5.36 

7.36 

10.34 

15 .20 

r ‘"•U 

0 

2 

3 

10 

5.59 

5.58 

6.99 

6.93 

9.03 

9.02 

12.12 

12.12 

17-15 

17. is 


0 

o 

3 

i0,79 

11.50 

13.04 

15.83 

20.76 


3 

10 

10,79 

11.50 

13.04 

15.83 

20.75 


0 

0 

8 

29.25 

39.79 

52.57 

72.08 

104.13 


10 

29.27 

39.80 

52.58 

72.21 

104.44 


0 

1 

8 

33.00 

41.71 

54.16 

73.44 

105.29 


10 

33.02 

41.70 

54.17 

73.57 

105.61 

c-c 

0 


8 

40.66 

41»e6 

59.04 

74.53 

108 .82 


2 

10 

40.64 

47.61 

59 .05 

77.68 

109. i4 




8 

53.17 

57 ,77 

67.31 

84.45 

114.74 


0 

3 

10 

53.15 

57.71 

67 . 34 

84.63 

115.07 





Table 

4, 
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V°r = » 

= 0.35 

rB' r 


0.3 Ne = 8 

B-C 

n 

m 









0.1 

0.2 

0.3 

0.4 

0.5 



0 

27.31 

34.61 

45.37 

62,00 

89,16 




27.30 

- 

45.20 


89.20 



1 

28.95 

36.11 

46.67 

63.13 

90,14 

C-C 

Q 

n 

28.40 


45.60 


90.20 



2 

36.65 

41.83 

51.18 

66.82 

93.22 




36.70 


51.00 

— 

93.30 



3 

51.25 

5 3.33 

60.09 

73.80 

98,80 




51'.2.f ■ . 


60.00 


99.00 



0 

17 .80 

22.71 

29.98 

41.26 

59.82 




17.80 

- 

29.90 


59.80 



1 

19.42 

24.27 

31.40 

42.56 

60.99 

<5— r* 

0 

3fr 

19.00 

- 

31.40 


61.00 


2 

26.74 

30.09 

36.25 

56*72 

64.64 




26.00 

- 

36.20 

- 

64,60 



3 

40.08 

41.32 

45 . 47 

54*38 

71.10 


— 

3t 

40,00 


45.40 


71.00 



0 

4.23 ' 

5.18 

6,66 

9.02 

13.02 



3f 

4.23 

~ 

6.66 

- 

13.00 



1 

3.48 

4.81 

6.55 

9,11 

13.29 

F~C 

0 


3.14 

- 

6.63 

- 

13.^ 



2 

5.62 

6.45 

7.95 

10.46 

14,70 




5.62 

— 

7,95 

- 

14.70 



3 

12.45 

12.61 

13.27 

14.96 

18,56 

— 


Jt 

12.40 

*• - 

13.30 


18 .50 


Table 5 
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Verification of the convergence of frequency parameter 

X ^ for orthotropic plates 

with linearly varying thickness 

^rp'^^r = °*35 = 0.3 

h^/h^ =1.1 


B-C n m Ne 






0.1 

0.2 

0.3 

0.4 

0.5 



0 

8 

3.28 

3.86 

5.20 

7.30 

10.81 



10 

3.27 

3.86 

5.20 

7.29 

10.81 



1 

8 

2.88 

3.97 

5.45 

7.65 

11.25 




10 

2.87 

3.96 

5.45 

7.65 

11.25 

C-F 

0 

2 

8 

4.14 

5.17 

6.68 

3.97 

12.69 




10 

4.14 

5.17 

6,68 

8.96 

12.69 



% 

8 

7.99 

8.51 

9.65 

11.72 

15.36 




10 

7.98 

8.51 

9.65 

11.71 

15.36 



0 

' 8 

24.03 

29.45 

38.90 

5 3.34 

77.07 



10 

23.98 

29.46 

28.87 

53.48 

78.12 


0 

■t 

8 

2 4.72 

30.83 

40.08 

54.36 

77.94 

G-C 

.u 

10 

24.67 

30.88 

40.05 

54.49 

79.00 



2 

8 

30.13 

35.27 

43.69 

57.38 

80.5 3 



<a 

10 

30.10 

35.27 

43.64 

57,54 

81.67 



o 

8 

39.35 

42.74 

'49.81 

62.50 

84.89 




10 

39.34 

42 .72 

49.75 

62.69 

86.16 





h./h^ = 1 

.2 





0 

d 

3.18 

z: 

3.39 

4.57 

6.40 

9.48 



10 

3.17 

3.39 

4.57 

6.40 

9,48 




S 

2.56 

3.48 

4.78 

6.71 

9.87 

C-F 


1 

r\ 

10 

8 

2.55 

3.64 

3.48 

4,55 

4.78 

5 ,86 

6.7i 

7,87 

9.87 

ii,i3 



2 

10 

3.63 

4.53 

5.86 

7.87 

11.13 




■ 8 

7.01 

7,47 

3.47 

10.28 

13,48 

. ... 


o 

10 

7.01 

7.47 

8.47 

10.28 

13.48 



0 ' 

■ 8 

22 . 12 

35.86 

34,15 

46.81 

67.64 



10 

22.05 

25.85 

34,13 

46.92 

67.58 



-i 

8 

21.83 

27.09 

35.18 

47.69 

68.40 



1 

10 

21.75 

27.09 

35 .17 

47 .80 

68 . 34 

c-c 



8 

26.47 

30.94 

38.34 

50.35 

70.68 


2 

10 

26.43 

30.95 

38.32 

50.47 

70.62 




8 

34,53 

37.51 

43,71 

54.85 

74.51 



3 

10 

34.54 

37.48 

43.70 

54.99 

74.45 


Table 6. 
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Frequency parameter X (I2p (l-v^ 
of a S~C plate 


D^/D 

Ne = 

3 D 

ris 

/D = 0* 
> r 

35 = 

Zf£:: 

0.3 



- ©■ r 

0.1 

0.2 

0.3 

0.4 

0.5 


0.8 

14.47 

21.66 

29,39 

40.73 

59.32 


0.9 

15.67 

21.85 

29.43 

40.76 

59.53 


1.0 

16.86 

22.04 

29.45 

40.32 

59.45 

0,4 

1.1 

18.04 

22.24 

29.50 

40.71 

59.71 


1.2 

19.21 

22.43 

29.50 

40.75 

59.14 


O.S 

14.90 

21.99 

29.65 

40.95 

59 .50 


0.9 

16.13 

22.18 

39.70 

40.98 

59.71 


1.0 

17.34 

22.38 

29.72 

41.04 

59.63 

0.7 

1.1 

13.55 

22,53 

29 .77 

41.02 

59.90 

.... 

1.2 

19 .75 

22.77 

29.77 

40.97 

59.32 


0..8 

15.31 

22.32 

29.91 

41,46 

59.69 


0.9 

16.56 

22.51 

29.96 

41.20 

59.90 

1.0 

1.0 

17.30 

22.71 

29.93 

41.26 

59.82 


1.1 

19 .04 

22.92 

30.03 

41.24 

60.09 


1.2 

20.26 

23.11 

30.03 

41.19 

59.51 


0.8 

15.71 

22.64 

30.17 

41.38 

59.38 


0.9 

16.98 

22.84 

30.22 

41.41 

60.09 

lo3 

1.6 

18.25 

23.04 

30.24 

41.48 

60.01 


1.1 

19.51 

23.25 

30.29 

41.46 

60.28 


1.2 

20.76 

23.44 

30.29 

41.41 

59.70 


O.S 

16.09 

22.96 

30.43 

41.60 

60.07 


0.9 

17 . 39 

2 3.16 

30.48 

41.63 

60.23 

1.6 

1.0 

18.63 

2 3.36 

30,50 

41.70 

60.20 

1.1 

19.97 

23.57 

30.55 

41.68 

60.47 



- 1.2 

21.25 

23.76 

30.55 

41.62 

59 .88 


Table 7 
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Frequency parameter X r2 :i2 p (i^v v )/Ce 
of a C~C plate " ^ ® ^ ■ 


Ne = 8 D^/D^=0.35 v^=0.3 


& r 

hi/hj 






0.1 

0.3 

0.3 

0.4 

0,5 


0.8 

23.53 

33.53 

44.82 

61.47 

88.54 


0.9 

24,98 

33.76 

44.88 

61.51 

88.62 

0.4 

1.0 

26.36 

33.93 

44.90 

61.62 

88.35 


1.1 

27.72 

34.21 

44.93 

61.29 

39.29 


1.2 

29.08 

34.43 

44.95 

61.95 

38.04 


0.8 

2 3.98 

33.85 

45.06 

61.66 

88.69 


0.9 

25.42 

34.07 

45.12 

61.70 

88.77 

0.7 

i.o 

26.34 

34.30 

45.14 

61.81 

89.01 


1.1 

28.24 

34.5 3 

65.17 

61.83 

89.45 


1.2 

29.63 

34.75 

45.18 

61.73 

88.19 


0.8 

24.40 

34,15 

45.30 

61.85 

83.85 

0.9 

25 .86 

34.38 

45.36 

61.89 

88.93 

1.0 1.0 

27.31 

34.61 

45.37 

62,00 

89.16 

1.1 

28.74 

34,85 

45.40 

62 02 

89.60 

1.2 

30.16 

35.07 

45.42 

61.92 

88.34 



0.8 

24.79 

34.46 

45.5 3 

62 .04 

89.00 


0.9 

26.28 

34.69 

45.54 

62 .08 

89 .08 

1.3 

1.0 

27 .76 

34.92 

45.61 

62.19 

89.32 


1.1 

29.23 

35 . 16 

45.64 

62 , 2 i 

89.76 


1.2 

30.68 

35.39 

45.66 

62.11 

89.49 


0.8 

25.17 

34.76 

45.76 

62.23 

89.16 

0,9 

26.70 

34.99 

45.82 

62.27 

89.24 

1.6 i.o 

28.21 

35.2 3 

45.84 

62.38 

89.48 

i.i 

29,70 

35.47 

45.87 

62.39 

89.92 

1.2 

31.19 

35.70 

45.89 

62.29 

89.65 


Table 8 
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Frequency parameter X rJ (12 P (1-v v_)/(e 

^ ^ r B r 2 

of a F-C plate 

Ne = 8 D„„/D_, = 0.35 V 


D /D =0.35 
na r 


9 


0.3 


D^/D^ h.As 


1.6 


0.8 

0.9 

1.6 

1.1 

1.2 


R1/R2 





0.1 

0.2 

I 

ro 

• 

0 

£).4 

0.5 



O.S 

\ 

2.68 

4.33 

5.98 

3.42 

12.47 



0.9 

3,03 

4.39 

5.99 

8.42 

12.47 

0. 

4 

1.0 

3.40 

4.45 

6.00 

3.42 

12.47 



1.1 

3.79 

4.52 

6.01 

8.42 

12,46 

- 


1.2 

4.18 

4.58 

6.03 

3.42 

12.47 



0.8 

3,08 

4.76 

6.322 

8.72 

12.74 



0.9 

3.45 

4.77 

6.33 

8.72 

12.75 

0. 

,7 

J.O 

3.84 

4.83 

6., 34 

8,72 

12.75 



1.1 

4.24 

4.89 

6.35 

8.72 

12 .75 



1.2 

4.66 

4.96 

6 . 36 

8.73 

12.75 



0,8 

3.43 

5.04 

■ 6.63 

9.01 

13.02 



0.9 

3.83 

5.11 

6.64 

9.02 

13.02 

1, 

.0 

1.0 

4.23 

5.18 

6.66 

9.02 

13.02 

-V ^ 


1.1 

4.66 

5.24 

6.67 

9.02 

13.02 



1.2 

5.09 

5.31 

6 .63 

9 .02 

13.02 



0.8 

3.75 

5.37 

6.94 

9.30 

13.28 



6.9 

4,16 

5.43 

6 .95 

9.30 

13.29 

1 

.3 

1.0 

1.1 

4.59 

5.04 

5.46 

5.57 

6.96 

6.97 

9.30 

9.30 

13.29 

13.28 



X JL 

1.2 

5.49 

5,64 

6.98 

9.30 

13.29 


4.04 

4.48 

4.93 

5.39 

5.86 


5.67 

7.22 

9.57 

5.74 

7.24 

9.57 

5.81 

7.25 

9.57 

5.88 

7.26 

9.57 

5,95 

7.27 

9.58 


13.55 

13.55 

13.55 

13.55 

13.55 


Table 9 
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edge free. 

IXae to the extra stiffness provided by the 
inplane stresses the value of a increases as expected. 

The variation of frequency parameter with different 
radii ratio is shown in Figs. 12 and 13 for different 
speeds: for the two boundari'' conditions, free-free and 
free-clamped respectively. The fundamental frequencv 
parameter is considered in both the cases. 

As in the case of free plates, the frequency 
decreases with radii ratio. The influence of rotation 
is more in the case of free-free disc than free-clanped 
disc. In Table 10 frequency parameters for free-clamped 
and free-free discs for various speeds and radii ratio 
are tabulated. In Table 11 convergence is verified by 
comparing the values obtained by 4 and 6 elements. The 
frequency parameters obtained using the itqilane stresses 
from P.E.M. and exact solutions are compared in Table 12. 

5,3d Rotating Ortho tropic Discs 

The variation of frequency parameter with speed 
is similar to that of isotropic plate, but for the change 
in magnitude. As in the case of isotropic plates it can 
be observed from Figs. 14, 15 and 16 that variation is 
more in the case of free-free discs than frae-claraped 
discs. The frequency parameters for various D^/D^ ratios 
and radii ratios are given in Table 13. 






Variation of the frequency parameter X R^CJ^h/D 

with angular velocity S for different radii ratio (R /R ) 

(inplane stresses from exact solu-fion) 


Ne = 5 


B-C 

Q 



Rj/Rj 





0.1 

0.2 

0.3 

0.4 

0.5 


100 

4.29 

5 .2 3 

6.71 

9.06 

13.06 


200 

4 ♦ 4 4 

5 . 38 

6.35 

9.19 

13.16 

■p-n 


300 

4,67 

5.62 

7.08 

9.40 

13.33 

m=0/ n=0 

400 

4.97 

5.94 

7.39 

9.68 

13.57 


500 

5.31 

6.31 

7.76 

10.02 

13.87 


100 

5.39 

5.24 

5.01 

4.71 

4.38 


200 

5,66 

5.52 

5.30 

5.03 

4.71 

P-P 

300 

6.07 

5.95 

5.76 

5.50 

5.17 


400 

6.61 

6,51 

6.34 

6.07 

5.71 

m=z / n=u 


500 

7,24 

7,16 

6.99 

6.69 

6.26 


Table 10 
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Verification of convergence of the frequency parameter 


^ r,2 

X 

( f h/D) for 

the rotatir^ di 

.3c (inplane stress 

from 

exact solution) 




V 

r 

=0,3 R^/R^ = 

: 0.5 

B-C 

R.P.M. 

m = 2,. 

n = 0 

Ne = 4 

Ne = 6 


100 

4,38 

4.38 


200 

4.71 

4.71 

P-F 

300 

547 

547 


400 

5.70 

5.71 


500 

6.26 

6.27 


100 

14.74 

14,74 


200 

14 86 

14-35 

P-C 

300 

15 .04 

15.04 


400 

15.30 

15.31 


500 

15.63 

15.63 


Table 11 


I 
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Comparison of the frequency parameter X (fh/D^)^/^ 
determined using inplane stresses from exact solutions 
and finite element methods. 



Q = 300, 

V = 0.3 , 

Ne = 14 

B-C 

RjAj 

F.E.M. 

Exact 


0.1 

4.74 

4.67 

F-C 

m=Oy 

0,3 

. '7*1^ 

7.08 

n=0 

0.5 

13.37 

13.33 


a.i 

6.30 

6,07 

F-i’ 

0.3 

5.94 

5.76 

m=2 

n=0 

0.5 

5.25 

5.1^ 


Table 12 




Frequency parameter 



Angular velocity (a) 


Fig. 15 Variation of frequency parameter of an 
orthotropic disc with speed 














Frequency parameter 
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Variation of frequency parameter X R^( h/D with 

angular velocity Q for 'different radii ratio and 

D^/D ratio 
fer XT 


= 180.00 D^/D^ =0.35 = 0.3 Ne = 12 


B-C R3_/R2 

© r 

Angular velocity Q 


100 

200 

300 

400 

500 


0.4 

3.46 

3.61 

3.86 

4.16 

4.52 


1.0 

3.30 

4.43 

4.75 

5.11 

5.51 


1.6 

4.99 

5.17 

5.44 

5.79 

6.19 



>fcNe=19 

0.4 

6.04 

6.16 

6.32 

6,57 

6.85 

F-C 

0.3 

1.0 

6.71 

6.87 

7.12 

7.47 

7.88 



1.6 

7.32 

7.51 

7.8 3 

8.24 

8.72 

(an.mj 


0.4 

12.49 

12.56 

12.70 

12.88 

13.10 

0^0 

0.5 

1.0 

13.06 

13. is 

13.38 

13.65 

13.99 


1.6 

13.60 

13.75 

14.00 

14.35 

14.77 




0.4 

4.02 

4.48 

5.15 

5.97 

6.88 


0.1 

i.o 

5.43 

5.79 

6.35 

7.06 

7.87 



1.6 

6.21 

6«55 

7.07 

7.73 

8.50 

F-F 

»Nee=19 

0.4 

3.47 

3.87 

ft. 48 

5.13 

3.39 


0.3 

i.o 

5.02 

5.35 

5.84 

6.47 

7,17 



1.6 

5.90 

6.26 

6.80 

7.4 

8.26 









0,2 


0.4 

2.90 

3.32 

3.88 

4.49 

5.09 

7 

0.5 

I.o" 

4.40 

4.77 

5.28 

5.87 

6.48 



1.6 

5.27 

5.67 

6.24 

6.91 

7.61 


Table 13 
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In Fig, 17 inplane stresses from F.E.M. and 
exact solntions are plotted. It is observed that conver 
gence is not so rapid as it was for vibration of plates. 
Hence higher number of elements (25) are required for 
better accuracy In case of rotating discs. However 
since stiffness due to inplane stresses forms a part of 
the overall stiffness, at low speeds, with lower number 
of d.ements ,(12) frequency parameter a is very close to 
the values obtained using the inplane stresses from 
exact solution. But it is obvious to go on increasing 
the number of elements till the values become almost 
stationary since convergence dop-cndn on otner factors 
boundary' conditions and radii ratio also. 








GHAPTrK ~ VT 


CONCLUSIONS 

6.1 ODNCLUSIONS 

Using semi analytical method vibration analysis 
of circular and annular plates for various boundary 
conditions has been done, both for isotropic and ortho- 
tropic plates. In the present analysis the problen of 
material singularity occurring in orthotropic plates is 
avoided by assuming a small hole at the centre and even 
then the answers are in good agreement with that of solid 
plates. The same element has also been used for the 
study of dynamic b^aviour of annular plates with varying 
thickness. In case of orthotropic plates it has been 
observed in all cases that the frequency parameter is 
mainly dependent upon the regidity modulus in radial 
direction compared to the material properties in other 
directions. It is because of this reason that when the 
frequency parameter is determined for a constant with 
varyirq frequency parameter does not vary much. Ihe 
conclusion to be drawn here is that, by providing suffi- 
cient reinforcement in radial direction (polar—orthotro— 
pic plates), frequency parameter as high as that for an 
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isotropic material can be obtained and these orthotropic 
plates are prefered due co their higher stiffness to 
weight ratio. 

In case of rotating discs the influence of D_ 
is more significant. So » for a rotating disc with cons- 
tant D^/when D^/D^ is increased/ it is found that freguency 
parameter increases more rapidly than that of a stationary 
plate. This is due to the obvious reason that inplane 
stresses increase as the regidity modulus increases. 

6,2 SCDPE FOR FUTURE WORK 

Since the semi analytical method is quite power- 
ful giving very good results (required for all engineer- 
ing purposes) further application of the method should 
be directed towards practical problons like rotating 
discs with periferal load (turbine disc with blades)and 
thermal stresses. 

Secondly in case of orthotropic plates with 
vairying thickness exact results arc not available. So 
it is necessary to determine the frequency parameter 
through experiments and conpare with the finite elatient 


method results. 
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APPENDIX-A 


BSAJ' SHAPE FUNCTIONS ARE AS FOLLOWS [23] 


N 
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N 
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1 
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4 0 ^ 
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APPENDIX-B 


B.l 

[r] matrix of 

EQN. (3.2) 

IS [ 11 ] 




a^(a-b) 

a^b 

b^( 3a-b) 

■w 2 

a b 



1 3 

72 


"”^T" 



6 ^a b 

-'a(a4-2b) 

" 6 a b 

-b( 2 a+b) 

1 1 

r= 

1 j 


~ 7 

.^2 



- 3(a+b) 

2 a+b 

3 (a+b) 

(a 4 - 2 b) 



~3 ■“ 



l 2 



2 

-1 

zl 

-1 



l3 

o' 

1 



J 

where 1 

= a - b 




B^? 

COEFFICIENTS OF 

OF EQN. ( 3 . 8 ) 


^11 

= (D^ + 4 D 

fcf 

re 

3 


^12 

= Dq (m - m ) 

^-2 



^13 

= (-2Vg 

/ 4 

+ Dq (m 

- 2r»'> - ^ “re 

m^) P_^ 

^14 

= (“SVgD^ + 

Dg (m* - 

3m^) - 3 

m^) P 

0 

^22 

= D_ (m“ - 1)^ 

k} 

^-1 




k 


23 


^24 


'33 


- (-2V^ ~ 1) + (m^ - 1) {m^ - 2)) 

= (-6Vq (ra^ - 1) + Dq {m*^ - 4m^ + 3)) P^ 


( 4 D, 


(m^ - 2) + (m - 2) + 4 D^m 




■)P- 
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= (12D^ ~ 6 V^ D^,(n^ - 2) - 2 Vq (m^ - 3) 

+ (m^ ~ 3) (m^ ~ 2) + 8 D^ m^) P 2 

^44 = - 12Vq D^ (m2 ~ 3) + D^ (m^ ~ 3 )^ 

+ 16D^q m 2 ) 

where 

^ 3 i 

P. = / h- r a.r 
b 

For linearly varying thichness plate i.e. 

h = oc + /3r 
these csDnstants are 

P __2 = (a 2 -h 2 )/ 2 ab + |32(a-b) + (a-b)/ab 

+ 3 a ^2 log (a/b) 

P 2 = (a-b)/ab + (a2-b2)/2 t log{aA>) 

+ 3 a ,02 (a-b) 

P ^ = a~ log (a/b) + (a2-*b2)/3 + 3a 0 (a-b) 

-1 

+ 3a 32 (a 2 -b 2)/2 

P - a2 (a-b) 4^^ (a^-b^)/4 + 30.^^ U^~h^)/2 
o 

+ 3a ^2 (a2-b3)/3 

P = a2 (a2-b2)/2 + 0^ (a^-b^)/5 + 3a2 ^ (a2-b^)/3 

+ 3a ^2 (a'^-b^)/4 
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Pj = (a’-b^)/3 + ^ (a^-b^)/4 

+ 5a (a^~b^)/5 

(a^~b‘^)/4 + jS^ (a’^~b*^)/7 + 3a^ ^ {a^~b^)/5 
+ 3a iS^ (a^~b^)/6 

h. a - hj ^ b, “ h. 

where a = —i ?; , ^ -= — — 

a - b a - b 

These relations can be used for uniform plate toy putting 
3 *s o in the above equations. 


B, 3 


ODEFFICIENTS OF 



OP EQN. 


(3.10) 


[m] 


(e) 


^1 ^2 ^3 


Sym, 


“3 <=4 


Q. 


where 

® i 

Q . = p / h r dr 

X 13 



For a plate with linearly varying thictoess i.e. 
h es a + 3r/ these coefficients are 


. ^ ia U^-hh/2 + 0 (a^-b3)/3) 

* f (a (a^-'bV3 ^ 

Q, « Y (a (a^-bV4 <aW)/5) 


■i 
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= ‘> (a (a^-b^)/5 + P (a^-b^)/6) 

Q- =5^ {a^-b^)/6 + 3 (a'^-b'^)/7) 

Qg = 9" (a^-lcP)/! + 3 (a^-b®)/8) 

= 9 (a (a^-b®)/8 + 3 (a^-b^)/9) 

These cx5eff icients can be used for uniform thichness 
plate by substituting 3 = o in the above equations. 
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APPE:IDIX-G 


C.l 


COEFFICIENTS OF OP EON. (4,2) 


= / [ 0^ (s, J [s ] + <^9 /r^ fs| [s]] 

P b '• 


In r <3iv 


For a uniform fhiclcness plate 




(e) 


where 

Q_. = 

Qo 

Q. 


Q 


2 
^3 

^4 

Qc 




Svm 


Q, 


Pj-Wj 


Q. 


®2+°3 

4P3-rt3 


Q, 


n 


3P3+Q3 


6P4+Q4 


9P-+QC 
D 3 


= h (Dj log (a/b) +03 (a~b)) 
= h (D. (a-b) + ^2 


(a^~b^)/3) 

(a4-b^)/4) 


h (Dj {a^-b ^)/2 + O 2 

h (D^ (a^-b^)/3 + ^2 

h (Dj (aW)/* + E2 (a®-b®)/5) 

h (Dj (a®-b®)/5 +^2 (a®-b®)/6) 


h 
h m' 
h m' 
h m' 


(D, (a®-b®)/6 + % (a''-b’)/V) and 


2 CE, i^^-bh/2 + % (a3-b3)/3),;.vj; 


2 (E, (a3-b3)/3 + ^ 

? <E, (aW)/4 + h 

1 


I V 4 ' i 
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h (E^ (a^-b^)/5 + (a^-b^)/6) 

h-m^ (E^ (a®-b®)/6 + (a'^-b'^)/?) 


where 


P'S, ' 


, Dj _ 


a - b 


<J~r ~ CTJ- 

1 2 ^1 ^2 
a~-~b'~ ■' “'2 ~ 'i[*'’3~b™ 


C.2 COEFFICIENTS OF [k^] OF, EQN. (4.9) 

k. = Cj- log (a/b) 

11 ^ 

k^ = ^*^5+^6^ 

J. 2 

k. = (2C. + C_) (a^-b^)/2 

"■13 ® 

k- = (3C^ + C^) (a^~b^)/3 

"■14 ' ^ 

\ ; ^ ^4 + c^) (a^~b^)/2 

22" 

ki = (2C^ + Cg + 3Cg) (a^-b^)/3 

2 3,', 

k. = (3C. + C_ + 4C^) (a'^~b‘*)/4 


4 ^4, 


K- = + ^5 4 (a ->b^)/4 


k^ = (6C^ + Cg 4- 5Cg) {a^-b''}/5 

3 4 

k. = (9C. + C- + 6C.) (a^-b®)/6 

i^4 4 o 6 


ki =k 

mn nm 
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where 


2T E h 
r & 



2'^ Eq h 

(1-Vr Vq) 



=4 ''e 





